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Many crucial biological processes are regulated by mechanical stimuli. Here, we report new findings that 
pico-Newton forces can drastically affect the stability of the site-specific DNA binding of a single 
transcription factor, the E. coli integration host factor (IHF), by stretching a short —150 nm DNA 
containing a single IHF binding site. Dynamic binding and unbinding of single IHF were recorded and 
analyzed for the force-dependent stability of the IHF-DNA complex. Our results demonstrate that the 
IHF-DNA interaction is fine tuned by force in different salt concentration and temperature over 
physiological ranges, indicating that, besides other physiological factors, force may play equally important 
role in transcription regulation. These findings have broad implications with regard to general 
mechanosensitivity of site-specific DNA bending proteins. 

Site-specific DNA binding by proteins is a crucial process for gene transcription and its regulation. For 
example, RNA polymerase binds to specific promoter sequences to initiate transcription, while many DNA- 
binding transcription factors also bind to these sequences in order to regulate the level of the transcription 1 . 
In eubacteria, e.g. E. coli, there are approximately 400 DNA-binding transcription factors 1,2 . Site-specific binding 
is also essential to many other DNA transactions, such as site-specific recombination and replication 3 5 . 
Therefore, detailed knowledge of DNA binding properties of site-specific proteins will provide important insights 
to their regulatory functions. 

It has been known that the activities of many site-specific DNA binding proteins are regulated by various 
environmental factors, such as temperature, salt osmolarity, pH value, as well as other cellular proteins. While the 
sensitivities of site-specific DNA binding to these factors have been extensively studied by various experimental 
approaches, the potential role of force has been much less investigated. 

Only very recently the mechanical force has been recognized as an important physiological factor that plays 
critical functions in diverse biological processes. In vivo, forces are produced by various cellular machineries. For 
example, in nucleus, individual DNA and RNA polymerases can exert up to —30 pN on DNA during actions 6,7 . 
Due to the ubiquitous presence of forces on chromosomal DNA, force may potentially influence site-specific 
DNA binding by proteins. In this work, we seek to probe the site-specific DNA-protein interaction on a single- 
DNA tether placed under force. As many DNA binding transcription factors cause DNA bending upon binding, 
the model system we chose is the E. coli integration host factor (IHF); a well-known DNA bending protein which 
binds to several specific sequences. 

IHF was originally discovered as an essential co-factor for site-specific recombination of phage X into E. coli 
genome 3 . It also regulates various gene transcriptions in E. coli" and S. typhimurium 9 . Various functions of IHF 
rely largely on its site-specific binding 1,10 . Among several specific sequences to which IHF binds with high 
affinity 4 , the —34 bp H' sequence in the bacteriophage X attachment site has the strongest binding affinity to 
IHF. Its dissociation constant is in the range of 2-20 nM 4 according to previous bulk biochemistry measurements, 
and <1 nM according to a recent stop flow measurement 11 . The IHF induced bending angle of H' was deter- 
mined to be —160° by X-ray crystallization 5 , >120° by atomic force microscopy (AFM) imaging 1214 . 

Besides site-specific binding, as an abundant nucleoid-associated protein (NAP) with an intracellular concen- 
tration of 12-55 uM 15 , IHF also associates with DNA non-specifically 16,17 . It is reported that in the presence of 
MgCl 2 in mM range, at high enough concentration, non-specific binding of IHF can mediate condensation of 
large DNA 16 , suggesting a role of non-specific binding in packaging and shaping the genome. 

Binding of IHF to specific sequences such as H' sequence has been investigated by electrophoretic mobility 
shift assays (EMSAs) and footprinting 4,18 , fluorescence resonance energy transfer (FRET) 19,20 , isothermal titration 
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calorimetry (ITC) 17 , stop-flow fluorimetry 11 , tethered particle 
motion 21 , and temperature jump 22 . IHF is able to bind to specific 
DNA sequences in both specific and non-specific modes. While 
non-specific binding may occur at high protein concentrations and 
low KC1 concentration, specific binding dominates at low protein 
concentration or high KC1 concentration 17 . In addition, several inter- 
mediate conformational states have been suggested by recent stop- 
flow fluorimetry 11 and time-resolved FRET 19 . 

However, in these studies DNA was not subject to force; therefore 
the potential role of force on specific IHF-H' interaction remains 
unclear. To address how force may regulate the interaction between 
IHF and H' sequence, here, we investigated the force-dependent 
stability of the conformational states of a single IHF-H' complex 
in different salt, temperature and IHF concentration conditions. In 
addition, we also determined the bending angle of H' induced by IHF 
in solution. 

Results 

Our study was based on DNA extension changes under tension 
induced by the docking of a single IHF protein to an H' sequence 
inserted to the middle of a 445 bp (—150 nm in contour length) 
DNA that does not contain any other high-affinity site based on 
the known consensus sequence for IHF (SM-Text "DNA con- 
structs") (Fig. la). The extension change was measured by magnetic 
tweezers with a spatial resolution of ~2 nm and a sampling rate of 
100 Hz 23 . To ensure the specificity of the binding, parallel experi- 
ments were performed on a control DNA of the same sequence but 
lacking the H' site (SM-Text "DNA constructs"). In E. coli, the 
MgCl 2 concentration is tightly regulated in a narrow range of 1- 
4 mM 24 . Therefore, in all our experiments of the specific interaction 
between IHF and H' DNA sequence, 2.5 mM MgCl 2 was included to 
mimic the physiological condition. 

Fig. lb shows representative time traces recorded in 10 nM IHF, 
50 mM KC1, 2.5 mM MgCl 2 , 10 mM Tris (pH 7.4), and 21 ± 1°C at 
force range from 0.5-1 pN. Clear two-state fluctuations between a 
shorter DNA extension (bent DNA) and a longer DNA extension 
(unbent DNA) were observed. The two-state fluctuations were fil- 
tered with fast Fourier transformation, followed by two-state digit- 
izing using a noise-beating step-finding algorithm similar to that 
developed by Cui et al 25 . (red data in Fig. lb) (see details in supple- 
mentary materials). Such two-state fluctuations were not seen before 
IHF was introduced (Fig. SI). It was also absent with the control 
DNA construct of the same sequence but lacking the H' site in the 
presence of IHF (Fig. S2). In addition, the force extension curves of a 
48,502 bp long X-DNA in the presence of 10 nM IHF were indistin- 
guishable from that of the naked X -DNA in 50 mM KC1 (Fig. S3a), 
which suggested that non-specific binding of IHF does not play an 
observable role in this condition (note that on lambda DNA there are 
four high affinity consensus sites 26 , but their contributions to the 
force-extension curves of the vast lambda DNA is negligible 27 ). 
Thus, we conclude that the two-state fluctuation was both IHF and 
H' sequence dependent. 

Fig. lc shows the normalized histogram of DNA extension dis- 
tributions, which are apparently bimodal with two extension species, 
fitted with double Gaussian distributions at different forces (Fig. lc). 
The distance between the two peaks calculated based on 40 such 
distributions is 17.34 ± 0.08 nm (mean ± s.d.) represented the 
extension difference of the two conformational states at the force 
range. Remarkably, a small force change of sub-pN switched the 
lower extension state from more probable to less probable, suggest- 
ing that the balance between the two states can be fine-tuned by 
adjusting force over a narrow range. 

The force-dependent probability of the bent DNA state, Pbenb can 
be calculated by the relative area of the two species in the bimodal 
extension distribution. This probability is related to the free energy 
difference between the two states through the Boltzmann distribution: 
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Figure 1 | Two conformational states of DNA-IHF complex, (a). A sketch 
of the principle of the measurement. The FT DNA bending by IHF binding 
reduces the DNA extension, which is detected by the resulting change of 
the height of the bead. (b). Dynamic fluctuations between two extensions 
in 10 nM IHF, 50 mM KC1, 2.5 mM MgCl 2 , 10 mM Tris (pH 7.4), and 
21°C at different forces (0.5-1 pN) indicated by different colors. The red 
line represents extension steps detected a noise-beating step-finding 
algorithm (SM-Text "Noise-beating step-finding algorithm"), (c). The 
probability densities of the two extension states, which were produced by 
the double Gaussian fitting of the normalized histogram of smoothed data 
at different forces (0.4-1 pN, represented by different colors) using a bin 
size of 1 nm. A statistics of 40 distributions gave a step size of 17.34 ± 0.08 
(mean ± s.d.) nm and R 2 = 0.96 ± 0.022 (mean ± s.d.). Note, only 13 
fitting lines from 5 DNA molecules were shown in figure for clarity. 



^bent — 



Exp[(AG-Azx/)/fc 5 r] 
l + Exp[(AG-Azx/)//c B T] 



(1) 



where AG is the Gibbs free energy difference between the unbent and 
bent DNA conformations in the absence of force, Az is the step size of 
two-state fluctuation, and/is the force applied to the DNA. While AG 
describes the stability of the interaction between IHF and H', it 
depends on the protein concentration C. The intrinsic stability of 
the interaction is described by the dissociation constant K^, which 
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is a concentration independent quantity and is related to AG through: 
K d = C X Exp( — AG/k B T). A lower value of K d indicates a stronger 
interaction. Substituting AG with K d and C, Eq. (1) becomes: 



Pbent — 



{C/Kd) xExp{-Azxf/k B T) 
1 + {C/K d ) x Exp( - Az xf/k B T) 



(2) 



This equation explicitly expresses the dependence of Pbent on both 
force and protein concentration, as well as the dissociation constant 
of the specific interaction between H' and IHF. In our experiments, C 
and /were given as experimental parameters and the resulting Pbent 
and Az were measured; therefore K d is the only free parameter that 
can be determined by fitting the experimental data to Eq. (2). As the 
fluctuation step size Az is narrowly distributed around its average, for 
simplicity Az was set to be a constant at its average value (17.34 nm) 
in our fitting. 

Fig. 2 shows the fitting at two different temperatures 21 ± 1°C 
(dark gray) and 31 ± 1°C (red) in 50 mM KC1, which yielded the 
average values of K d associated with fitting standard errors of 0.44 ± 
0.06 nM and 1.54 ± 0.08 nM, respectively, using data from multiple 
(£3) independent experiments under respective temperatures. The 
results showed that the bent DNA state in an IHF-H' complex 
became less stable as temperature increased. 

At a higher salt concentration of 200 mM KC1 containing 2.5 mM 
MgCl 2 , two-state fluctuation was not observed at 21°C, indicating 
that the unbent DNA state dominated in the force range in the 
presence of 10 nM IHF (Fig. S4). We reasoned that this could be 
due to reduced binding affinity at higher salt concentration; there- 
fore, we investigated IHF binding at higher IHF concentrations of 
100, 500 and 1000 nM, where two-state fluctuations were observed 
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Figure 2 | The effects of temperature, KC1 concentration and IHF 
concentration on H' DNA bending. The bending probability as a function 
of force in 2.5 mM MgCl 2 , 10 mM Tris (pH 7.4), at 10 nM IHF, 50 mM 
KC1, 21° (dark gray symbols), 10 nM IHF, 50 mM KC1, 31°C (red 
symbols), 100 nM IHF, 200 mM KC1, 21°C (blue symbols), 500 nM IHF, 
200 mM KC1, 21°C (orange symbols), and 1000 nM IHF, 200 mM KC1, 
21 °C (wine symbols). Data for each solution condition were obtained from 
multiple (S3) independent DNA molecules. Error bars for each data 
points (symbols) were standard deviations from multiple (S3) repeating 
measurements for the same DNA molecules. The bending probabilities 
were calculated by the relative area of the two species in the bimodal 
extension distribution. Under each solution condition, data obtained were 
fitted by the two-state model (Eq.(2)) to obtain and the standard error 
of K d (fitting error), which are indicated in figure panels by corresponding 
colors. The goodness of the fitting (R 2 ) are 0.84, 0.98, 0.94, 0.93, 0.99, 
respectively. The purple dot line is the theoretical calculation of the 
bending probability in 10 nM IHF, 200 mM KC1, 2.5 mM MgCl 2 , based 
on the two-state model with an averaged = 28.5 nM. 



(Fig. S4). Fig. 2 shows the resulting force dependent bending prob- 
abilities (symbols), the corresponding two-state theoretical fitting 
based on Eq. (2) (lines), indicated by blue, orange and wine colors, 
respectively. K d estimated from the three different IHF concentra- 
tions at 200 mM KC1 were overall consistent with each other, with a 
value of 28.5 ± 3.9 nM (mean ± s.d.). Based on the estimated at 
200 mM KC1, the bending probability at 10 nM IHF in the force 
range of 0.5-1 pN was predicted to be lower than 5%, which 
explained why we did not observe two-state fluctuations in 10 nM 
IHF at 200 mM KC1. 

Here we note that the salt dependence of IHF binding to DNA is 
rather complex, related to release of counterions when IHF initially 
weakly binds to DNA and further release of counterions which is 
countered by simultaneous release of water molecules in the follow- 
ing tightly bending process. As a result, the first step is sensitive to 
KC1 concentration while the second is less sensitive 28,29 . Increasing 
KC1 concentration will reduce the first step of binding; therefore 
higher IHF concentration is needed to observe the two-state IHF 
dependent fluctuation in our experiments. 

The lower extension state can be explained by bending of the H' 
DNA sequence associated with an IHF. However, the longer exten- 
sion state may have two alternative explanations: a DNA without an 
IHF bound at H', or an IHF-H' complex in a non-bent, intermediate 
state. The values of K d of ~1 nM for 50 mM KC1, 2.5 mM MgCl 2 , 
and that of ~28 nM for 200 mM KC1, 2.5 mM MgCl 2 determined in 
our experiments agree with values determined in previous bulk 
experiments 411 , suggesting that the two-state fluctuations in our 
experiments are likely dominated by dissociation and association 
of IHF to the H' sequence of DNA. 

The specific structure of IHF-H' complex was solved by X-ray 
crystallization, which indicated that the 34 bp H' DNA is bent over 
the surface of an IHF hetero-dimer by an angle ~160° 5 . The H' 
bending angle by IHF was also estimated by AFM imaging to be 
>120 ol2 ~ 14 . Both of these methods determined the angles based on 
static conformations of the IHF-H' complex. The DNA bending 
angle can also be estimated from our dynamic two-state extension 
fluctuation data under constant forces in solution. DNA can be 
described as a semi-flexible worm-like-chain (WLC) polymer model 
with a DNA bending persistence length of A ~ 50 nm 30,31 . In the 
WLC model, a DNA of contour length of L can be modeled as a chain 
of N segments, each with a segment length of b <C A. The bending 
energy cost in units of k B T of one DNA segment is described by: 



E,m +1 ) 



IA 

2b 



(ft 



A s2 
-H+l) > 



where i denotes the ! tn vertex connect- 



ing the z' tn and (i + 1 ) th segments, and f; is the tangent vector of the ; th 
segment. The total bending energy is the sum of all the bending 

energies carried by each vertex: E= £;(£), r) + i). 

In our experiments, the 445 bp DNA contains one specific H' site 
near the middle of the DNA that is subject to binding and unbinding 
of IHF. To estimate the level of overall extension reduction of the 
DNA tether under constant forces when an IHF is associated to the 
H' site, we discretized our DNA by a segment length of b = 1 nm. 
The H' site associated with an IHF was modeled as a point-like 
kinked site placed at the middle of the DNA, which has a preferred 
bending angle described by a parameter y: 0 = cos _1 y. The bending 
energy of this site is therefore modified to be: 



-Ekink (»,*') 



-y) 



where the i and f are the tangent vectors of the 



two successive segments involved in the kink; the dimensionless 
parameter a describes the deformability of the kinked site, which 
was chosen to be a large number to ensure a rigid kink (a = 
50 nm was chosen in our calculation for Figure 3, but note the 
extension reduction versus force profiles are insensitive to the value 
of a as shown in Figure S6). The rest of DNA remained in the naked 
DNA state and follow WLC model. 
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The force-extension curves of the DNA tether with various values 
of preferred bending angle 0 at the kinked H' site were calculated 
using a previously developed transfer-matrix method 27,3233 . The 
curves shift downward as 0 increases (Fig. S6a), leading to an exten- 
sion reduction Az compared to DNA without the kink defect. Fig. 3 
shows Az versus force at several DNA bending angles in the range of 
90-180°, converted from y in the range of [— 1,0] by 0 = cos _1 y. In 
our experiment, the force range where the two -state fluctuations 
were observed was around 0.5-1 pN, which was associated with 
extension fluctuation step sizes (i.e, Az) of — 17.34 nm. By compar- 
ison with the theoretical prediction in Fig. 3, these values corre- 
sponded to predictions with a preferred DNA bending angle of 
140-180°. The bending angle 0 was also independently estimated 
by an approximate analytical formula derived by Kulic et al 34 .: 

Az = 4(1 — cos | — ) 1 4 / — , which is valid for forces 



k R T 



f 



^0.08 pN where the entropic DNA conformational fluc- 



tuation is suppressed. This analytical formula also gave similar esti- 
mated bending angle to that observed in our experimental force 
range (Fig. S6). Overall, the bending angle range estimated from 
our dynamic fluctuations is consistent with the large H' DNA bend- 
ing observed from previous X-ray crystallization and AFM imaging 
experiments 5,12 " 14 . 

Discussion 

In the present study, we showed that specific DNA-protein interac- 
tions can be probed by stretching short DNA molecules. Previous 
single-DNA stretching studies of DNA-protein interactions were 
usually restricted to DNA length of a few microns 35 , while some 
recent studies have pushed the length limit to —700 nm 
(—2000 bp) 36,37 . It was shown that the protein density on DNA must 
be high enough in order to see a detectable effect. As an example, for 
DNA bending by 90°, the spacing between two adjacent bends needs 
to be below 100 nm in order to see an obvious shift in force-extension 
curves 27 . The low signal-to-noise ratio resulted from large conforma- 
tional fluctuation of long DNA has made it difficult to detect protein 
binding at a low protein density on DNA 27 . As a result, those previous 
studies were mainly restricted to non-specific DNA-protein interac- 
tions at rather high protein concentrations to ensure sufficient bind- 
ing density for detection, including studies on bacterial IHF 16,38 , 




0.1 1 

Force (pN) 

Figure 3 | Theoretical prediction of the effect of bending angle and force 
on the extension of short DNA. The extension reduction (Az) of a DNA 
with a contour length of 1 50 nm induced by a kink placed at the middle as a 
function of force and bending angles 8 = 90° (gray dotted line), 143° (gray 
solid line), 162° (black dotted line), 180° (black solid line). The shadow 
area represents the rough force range and Az range measured in 
experiments. 



H-NS 39 - 41 , HU 42,43 , archaeal Alba, Cren7 and Sul7 44,45 , and eukaryotic 
HMGB 46 . 

Compared to previous biochemistry methods to probe specific 
DNA-protein interactions 4,11,17,18,22 , stretching short DNA has unique 
advantages to directly probe real time distortion of DNA by binding 
of a single protein. Compared to fluorescence based single-molecule 
method such as FRET 19 , single-DNA stretching does not require 
fluorescence labeling and is able to monitor protein-DNA inter- 
action over much longer time scales and distance range. Compared 
to the tethered particle motion that also probes the dynamics of DNA 
conformational changes by protein binding 21 , single-DNA stretching 
provides additional information of the potential regulatory role of 
force applied to DNA. Thus, our results demonstrate that stretching 
short DNA provides a powerful experimental approach for studies of 
specific DNA-protein interactions under force, which is compli- 
mentary to previously developed biochemistry and biophysical 
methods. 

The intracellular concentration of IHF (12-55 uM 15 ) is much 
higher than the concentration needed for site-specific functions; 
therefore non-specific binding of IHF also occurs. When IHF binds 
to DNA non-specifically, it occupies less number of DNA base pairs 
(—10 bp) 17 , and causes a much less DNA bending angle (<30°) at 
the binding sites 16 . This is in sharp contrast to specific binding of IHF 
to the H' sequence that occupies 34 bp DNA and causes a bending 
angle of >120° 5,12 - 14 . 

Despite that in vivo the complicated protein-DNA interactions 
may significantly differ from in vitro studies in kinetics and stability, 
our results of the ultra sensitivity of the IHF-H' to force provided 
important insights to the potential physiological role of force in 
regulating the activities of transcription factors. The nucleoid is 
known attached to the cell wall 47 , which may lead to building up 
tension in DNA due to DNA compaction by nucleoid associated 
proteins. Forces in picoNewton range can be expected based on 
A: B T interaction energy between proteins and DNA with nm scale 
of interaction distance. Force on DNA can also be produced during 
transcription, as RNA polymerases are known able to generate >20 
pN forces 6 . Therefore, our finding that IHF-H' interaction is sens- 
itive to sub-picoNewton forces is likely physiologically relevant. 

We showed that the binding and unbinding of IHF to its specific 
site can sense small forces < 1 pN, and can be fine tuned by sub pN 
force changes. At zero force, Pbent is larger than 95% at 50 mM KC1, 
2.5 mM MgCl 2 and 21°C (Fig. 2), consistent with the previous 
reported stable IHF-H' complex in the absence of force 4,11,12,21 . 
Such ultra force sensitivity of the IHF-H' interaction relies on a large 
dimensional change of DNA by protein-induced DNA bending. The 
— 17 nm step size resulted from large DNA conformational change 
observed in our experiments implies that a small force around 
k B T 

~ 0.2 pN can perturb the IHF-H' interaction, which is close 

17 nm F v 

to the force range observed in our experiments. As the functions of 
many site specific transcription factors rely on sharp DNA bend- 
ing 48,49 or DNA looping 50,51 that cause large conformational change 
of DNA, our results suggest that force on DNA may sensitively 
regulate the functions of a broad class of site specific DNA binding 
transcription factors. 

These results also demonstrated that the level of DNA distortion at 
the specific site bound with a protein can be estimated from the 
dynamic extension fluctuation between protein-bound and pro- 
tein-unbound states based on DNA polymer physics. In the case of 
IHF-H' interaction investigated in our work, the difference of the 
DNA bending angle (140-180°) estimated from our data and that 
determined in X-ray crystallization (— 160°) 5 is within 15%. A unique 
advantage of the approach of short DNA stretching to determine 
protein induced DNA distortion is that it is not limited to any specific 
solution conditions such as temperature, salt, pH, and protein 
concentrations. 
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Previously, IHF binding to a short DNA tether of 76 bp containing 
the H' sequence was studied by monitoring the displacement of a 
bead attached to one end of DNA from a cover glass in the presence 
of 100 nM IHF in lxPBS 52 . In that experiment, longitude fluctuation 
of the bead was recorded and used to re-construct the potential. 
Force was calculated by the gradient of the potential, which depended 
on the distance between the bead and the cover glass surface and 
therefore was not constant during DNA extension fluctuation. An 
~8 nm extension change between the bent and unbent DNA con- 
formation states was reported 52 , which was ~ Vi of that in our obser- 
vation. We propose that this difference might be caused by different 
level of DNA wrapping on IHF. Most of DNA bending by IHF 
heterodimer mainly occurs at two large kinked steps where the minor 
groove of DNA at each location is intercalated by one proline residue 
at the tip of the |3-ribbon arm from each IHF subunit 5 . In the work by 
Dixit et al, tension on DNA increased rapidly during the wrapping 
process 52 , which might prevent complete wrapping of DNA on IHF, 
i.e., possibly only one location is intercalated by a proline residue 
from one of the IHF subunits. This would lead to a smaller bending 
angle and explain why they observed a smaller extension difference 
between naked and wrapped DNA than we observed. In contrast, full 
DNA wrapping might occur in our experiments, as during wrapping 
a constant force is maintained. Our current study has further 
improved from that pioneering study by applying accurately con- 
trolled constant forces to DNA, which allowed us to extract more 
comprehensive information of the interaction between IHF and H' 
DNA, such as the dissociation constant and the bending angle of 
the IHF-H' complex. 

Several possible intermediate conformational states have been 
suggested previously 1119 . Using a noise-beating algorithm, some 
minority transitions associated with <10 nm steps were detected, 
which is consistent with the presence of such intermediate states (See 
Figure S7). Future work will be directed to improve the spatial reso- 
lution of our instrument as well as combine with single-molecule 
FRET technology to better resolve those possible intermediate states 
and their stabilities. 

Methods 

An improved magnetic tweezers with a spatial resolution of ~2 nm and a sampling 
rate of 100 Hz which is capable of stretching both short DNA and protein 23,53-55 was 
used in the experiments. Details of DNA construct and protein, magnetic tweezers 
measurement, temperature control, data analysis methods and supplementary figures 
can be found in the Supplementary information. 
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